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Abstract. Measurements of neutrino- nucleon and antineutrino-nucleon differential 
cross sections using the CCFR neutrino detector at Fermilab have been used to extract 
preliminary values of R = ol/o't in the kinematic region 0.01 < x < 0.6, and 
4 < Q 2 < 300 GeV 2 . The new data provide the first measurements of R in the 
x < 0.1 region. The x and Q 2 dependence of R is compared with a QCD based fit to 
previous data. The QCD fit, which provides an estimate of R in small x region where 
R has not been previously measured, is in good agreement with the new CCFR data. 
(To be published in Jou. Phys. G) 



1. Introduction 



The ratio R = ctl/ct of the longitudinal and transverse absorption cross sections in deep 
inelastic lepton-nucleon scattering provides information about the transverse momentum 
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and spin of the nucleon constituents. Within the theory of quantum chromodynamics 
(QCD), the nucleon constituents are spin 1/2 quarks and spin 1 gluons. In leading 
order QCD, R = 0, since the quarks have no transverse momentum. In the next to 
leading order formalism (NLO), to first order in a s , R is non-zero because of transverse 
momentum associated with gluon emission PJ. The NLO QCD prediction is given by an 
integral over the quark and gluon distribution and is proportional to a s . A measurement 
of R is a test of perturbative QCD at large x, and a clean probe of the gluon density at 
small x where the quark contribution is small. 

Poor knowledge of R , especially at small x, results in uncertainties in the structure 
functions extracted from deep inelastic scattering cross sections. The most precise 
previous measurements of R are from SLAC electron scattering experiments ||, and 
are at high x and low Q 2 , where non-perturbative contributions (e.g. target mass and 
higher twist effect) are important. Therefore, extraction of structure functions at large 
Q 2 or small x must rely on good measurements of R, or on reasonable estimates of R 
in the kinematic region where R has not been measured. 

We report here on a new extraction of R from measurements of neutrino and 
antineutrino differential cross sections in the CCFR neutrino detector. We compare 
the data (particularly at small x) to the prediction of a QCD based fit to R by Bodek, 
Rock and Yang || . Extrapolation of such a QCD based fit to smaller values of x can be 
also used to obtain estimates of R in the HERA region where R has not been measured. 

2. R calculation within QCD 

In this section we provide a brief description of the QCD based model of Bodek, Rock 
and Yang (BRY). They have used a calculation @ of the QCD contribution to R to 
order a 2 (NNLO). The contribution of massive heavy quark (charm) effects to R has 
been included within the photon gluon fusion process ||. These higher order QCD 
corrections and heavy charm quark effects are important at small x. In addition, the 
SLAC data on R indicate that there are additional non-perturbative 1/Q 2 and l/Q 4 
contributions at low Q 2 @, [| . One of these contributions, target mass effects, dominates 
at low Q 2 and large x. A comparison of the calculation (including both the heavy 
charm contribution to order a 2 and the Georgi-Politzer target mass corrections 0) 
to the SLAC data indicates that additional small higher twist contribution at low Q 2 
are required ||. A parameterization of the higher twist effects at low Q 2 should be 
constrained such that R approaches zero at the photoproduction limit (Q 2 = 0). BRY 
find that a simple empirical parameterization of the higher twist contribution to R (of 
the form A[(Q 2 - B)/(Q A + G)\) fits the SLAC data very well. The best BRY fit to the 
all previous data [0, |7| (for Q 2 > 0.4 GeV 2 ), using the GRV94 parton distribution 
for the QCD contribution, yields A = 0.37, B = 0.82, and C = 0.54 (x 2 /dof is 21.7/19). 
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Figure 1 shows some of the previous data and the model predictions calculated using 
various parton distributions. All curves include QCD to order a 2 s , the Georgi-Politzer 
target mass correction and the empirical small higher twist contribution. The curves 
labeled with (h) also include the effects of the charm quark mass. Also shown for 
comparison are the extrapolation of a previous empirical fit (Rworld) for R used by 
the SLAC, CCFR, and NMC collaborations. The QCD model predictions with various 
parton distributions || [| all yield similar results for R at large x (> 0.1). At small x, 
the predictions of the model using either the MRSD'_ and the MRSA distributions are 
very different because of the different gluon distributions. Note that target mass effects 
at small x are negligible, and the small empirical higher twist contribution is assumed to 
be independent of x. The model calculated with the GRV94 parton distributions (solid 
curve, figure 1) is probably the best available estimate for extracting F 2 in regions where 
R has not been measured. We proceed to test the models shown in figure 1 against the 
new CCFR data in the section below. 
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Figure 1. R at fixed x vs Q 2 : A comparison of the BRY prediction of R calculated 
using various parton distributions (including target mass and empirical higher twist) 
to data. A massive heavy (charm) quark contribution is included in the MRSA(h) 
and GRV94(h) prediction. Also shown is Rworld, which is a previous empirical 
parameterization of R data in the SLAC region. 



3. A Preliminary CCFR measurement of R 

The data sample includes two CCFR data runs (E744 and E770) collected using the 
Fermilab Tevatron Quad- Triplet neutrino beam. The wide-band beam is composed of 
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and with energies up to 600 GeV. The CCFR neutrino detector ]ID| consists of 



an unmagnetized steel-scintillator target calorimeter instrumented with drift chambers. 
The hadron energy resolution is AE/E = 0.85 /\^E '(GeV). The target is followed by a 
solid iron toroidal magnet muon spectrometer which measures muon momentum with 
a resolution Ap/p = 0.11. The relative flux at different energies is obtained from the 
events with low hadron energy, v < 20 GeV, and is normalized so that the neutrino 
total cross section equals the world average <j vN / E = (0.674 ±0.014) x 10~ 38 cm 2 /GeV, 
and a VN /a vN = 0.504 ± 0.005. 
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Figure 2. Preliminary differential cross section (statistical error only) for neutrino [a] 
and antineutrino[b] at E v — 150 GeV: Preliminary CCFR data (solid circles) are 
compared to the QCD LO prediction and previous CDHSW data (open squares). 



Preliminary values of R have been extracted from the differential cross sections in 
neutrino charged current scattering. The differential cross sections are written as 
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Edxdy 7T [ " V '\ a ' 2(1 + i2) 
where y = v / E u is the inelasticity, and R = a^/aT = F 2 (l + Q 2 /V 2 ) /2xFi — 1. In the 
quark-parton model, 2xF^ v = q + q and xF^ v = (q — q) ±2(s — c) on an isoscalar target. 

The data sample consists of 1,280,000 and 270,000 events after fiducial and 
kinematic cuts (P M > 15 GeV, ^ < .150, v > 10 GeV, Q 2 > 1 GeV 2 , and E v > 30 
GeV) . Dimuon events are removed because of the ambiguous identification of the leading 
muon for high y events. The data are corrected by Monte Carlo for cuts, acceptance 
and resolution smearing. The differential cross sections are determined in bins of x, 
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y and E v (0.01 < x < 0.60, 0.05 < y < 0.95, and 30 < E v < 360 GeV). Figure 2 
shows the differential cross sections at E v = 150 GeV. The differential cross sections 
agree with a LO QCD model which predicts a quadratic y dependence at small x for 
neutrino and anti-neutrino and a flat y distribution at high x for the neutrino cross 
section. A comparison of CCFR and CDHSW differential cross sections J7| indicates 
good agreement at high x, and poor agreement at small x, as shown in Figure 2. 

The values of R are extracted from linear fits to F vs e at fixed x and Q 2 bins 
through the following relation, 



F(x,Q 2 



Tr(l-e) 
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2xFi(x,Q 2 ) l + eR(x,Q 2 ) (2) 
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Here e ~ 2(1 — y)/(l + (1 — y) 2 ) is the polarization of 
virtual W boson. Corrections for electroweak radiative effects (Bardin), the W boson 
propagator, and isoscaler target (the 6.85% excess of neutron over proton in iron) are 
applied. The correction for AxF 3 = xF£ — xFg (which is important at small x) requires 
a knowledge of the strange sea and charm production mass effects (e.g. slow rescaling). 



The strange sea distribution, extracted to NLO from CCFR dimuon data [jTIJ has been 
used. The fits for R at each x and Q 2 have reasonable \ 2 ■ Figure 3 shows one of the 
linear fits. The values of R with statistical errors at fixed x vs Q 2 are shown in Figure 4. 
The new data, which extends a decade lower in x, are in agreement with the QCD based 
model. The values of R are not sensitive to the normalization of the absolute flux. The 
dominant uncertainty at very small x and low Q 2 is from the non-zero value of AxF$ 
because of uncertainties in the slow rescaling formalism and the level of the strange sea. 
Therefore, the values of R in the kinematic region x < 0.1 and Q 2 < 4 GeV 2 are not 
shown at present. 
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Figure 3. F(x,Q 2 ,e) vs e fit at x = 0.045 and Q 2 = 5.01 GeV 2 : The errors on F are 
only statistical errors. 



Better modelling of AxF 3 at low Q 2 and small x, and more precise data from the 
1996-97 run of NuTeV (FNAL E815) experiment with sign selected neutrino beams are 
expected to yield better determinations of R at small x. 



6 



PS 



PS 



PS 



PS 



0.75 
0.50 
0.25 
0.00 
-0.25 

0.75 
0.50 
0.25 
0.00 
-0.25 

0.75 
0.50 
0.25 
0.00 
-0.25 

0.75 
0.50 
0.25 
0.00 
-0.25 
-0.50 








0.5 1 


5.0 


10.0 

Q 2 


50.0100.0 


• 


CCFR 




— MRSD'_ 


O 


SLAC 




- MRSA(h) 


□ 


CDHSW 




* 


BCDMS 




GRV94(h) 


H 


EMC 




Rworld 



0.5 1.0 5.0 10.0 



Q 

Figure 4. R at fixed x vs Q 2 : Preliminary CCFR data are compared with other data 
and with the Bodek, Rock and Yang QCD based fit with various parton distributions. 
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